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Abstract. The well known scaling of the Edwards- Wilkinson equation is essentially 
determined by dimensional analysis. Once a drift term is added, more sophisticated 
reasoning is required, which initially suggests that the drift term dominates over the 
diffusion. However, the diffusion term is dangerously irrelevant and the resulting 
scaling in fact non-trivial. In the present article we compare the resulting scaling 
of the Edwards- Wilkinson equation with drift and Neumann boundary conditions to 
the published case with Dirichlet boundary conditions. 
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1. Introduction 

The Edwards- Wilkinson equation [TJ is very well understood. It is a stochastic equation 
of motion describing the most basic surface evolution of a growth model, consisting 
merely of diffusion and random particle deposition. Given its minimal parameterisation 
in terms of a diffusion constant D, a noise amplitude T as well as the system size L and 
time t, Equation ([T|) and ()2]), the scaling of its roughness in L and t is determined by 
dimensional analysis, because the only dimensionless quantity is Dt/L 2 . This changes, 
as soon as another interaction term is added whose coupling allows an additional 
dimensionless quantity to be formed. In principle, the scaling is then not determined 
by simple dimensional arguments. In this situation, it is common to invoke physical 
arguments to show that one of the terms is (infrared) irrelevant and thus can be dropped 
[2j [3j H]. There is no mathematically rigorous argument underpinning this procedure 
and it therefore can produce incorrect results. This has been analysed in the past (5] for 
the addition of a drift (convection) term to the Edwards- Wilkinson equation, which first 
seems to render the diffusion irrelevant, when in fact it becomes dangerously irrelevant 
[6] and thus should not be dropped from the analysis. 

It has been noted (5] that the role of the drift term depends on the boundary 
condition. It is trivial to show that in the presence of periodic boundaries the drift does 
not enter the roughness (in the definition below) at all, but changes the scaling when 
Dirichlet boundary conditions are applied. One might now wonder to what extend the 
result is universal, i.e. what happens, for example, if Neumann boundary conditions are 
applied. The calculations are far more involved than for the Dirichlet case and presented 
in the following. 

It is not surprising if boundary conditions change some universal quantities at the 
critical point [7] , certainly not if there is some net transport across the system as in the 
present case. Given that the system is apparently sensitive to the choice of boundary 
condition, the finding below might come as a surprise; Neumann and Dirichlet boundary 
conditions produce, in fact, the same leading order behaviour. 

Beyond its role of probing scaling and universality in the Edwards- Wilkinson 
equation, the drift term has a simple physical motivation. It could, for example, be 
caused by a gravitational or an electric field, present during epitaxial growth. The 
choice of the form of the boundary condition is a bit more arbitrary. A Dirichlet 
boundary condition corresponds to particles falling off the edge of the substrate. At 
the same time, particles are to be generated whenever a hole forms at the boundary, 
which makes the boundary condition somewhat unphysical. As discussed below, the 
Neumann boundary is an improvement in this respect, but still does not correspond 
an easily implemented physical situation, for example by allowing no net flux across 
the boundary. The motivation for the choice of the Neumann boundary condition in 
the following is thus the question to what extent the choice of the boundary condition 
changes the universal behaviour of the system compared to the Dirichlet case. 
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1.1. Preliminaries 

The features of the original Edwards- Wilkinson (EW) equation [I] 

^(x,t) = J DV 2 0(x,t) + r / (x,t) (1) 

have been reviewed a number of times [2, El EJ [El [H]. The field 0(x, t) describes 
the height of an interface over a substrate of dimension d and linear extension L at 
position x and time t. The noise i]('x,t) is Gaussian and white, characterised by the 
usual correlator 

(r / (x,t)r / (x',t')) = r 2 5(x-x')5(t-t / ) (2) 

with amplitude T 2 and (rj) = 0. The operation (•) describes the averaging over the 
noise, such that, for example, <9 t (</>(x, t)) = DV 2 (<^(x, t)}. In the context of growth 
phenomena, the observable of choice is often the roughness, which is the ensemble 
average of the spatial variance of the field 0(x, t) 

w 2 (L,t) = (Jt^tJ -Ji^f) (3) 



d d x ■ (4) 



where • describes the spatial average 
L d Jl* 

over the entire substrate of volume L d . The roughness is closely related to the correlation 
function C(x, x',t) = ^(0(x, t) — 0(x',t)) 2 ^>, as 

w2{L > t] = \h L ddx 'j? J L / dxC ^^ • ( 5 ) 

The scaling of the roughness in time is characterised by the universal dynamical 
exponent z and the finite size scaling in L by the universal roughness exponent a, which 
is summarised in the Family- Vicsek scaling hypothesis [12] 



w 2 (L,t) = aL 2 «g (A) (6) 

where a and b are two non- universal metric factors |7J. The (universal) dimensionless 
scaling function G(x) mediates between two regimes, as lim^oo Q (x) = Goo with 
< Goo < oo and G(x) ~ x 213 for x — > with the roughening exponent (3 obeying 
zfi = a. Taking the limit t — > oo first, the roughness thus scales like w 2 (L,t) oc L 2a in 
increasing L; taking the thermodynamic limit L — > oo first, it scales like w 2 (L,t) oc t 2/3 
with increasing t. Subleading terms in L and t respectively are expected to be suppressed 
with increasing t and L respectively. 

Given that the Edwards- Wilkinson equation in the form ([1]) is fully parameterised 
by L,t and D,T 2 , the only possible scaling is 

w 2 (L^D,T 2 ) = ^L 2 - d G^) (7) 

so that a = (2 — d)/2 and z — 2, unless w 2 does not exist, diverges or vanishes. In fact, 
at and above the upper critical dimension d c = 2, the roughness is controlled by a lower 
cutoff or lattice spacing a, and diverges for d > d c like a 2 ~ d as a — > [S]. 
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In the following, the Edwards- Wilkinson equation is studied in d = 1 dimensions, 
where 

a = 1/2 (3 = 1/4 z = 2 (standard EW) . (8) 

Provided that no other scales are present, the result (J7|) from dimensional analysis 
determines the scaling. Boundary conditions generally affect the scaling function Q(x), 
and the metric factors a and b, fl6]), but not the exponents. Some types of interactions 
can be added to the Edwards- Wilkinson equation without affecting its scaling behaviour, 
which hints at the universality alluded to earlier. This can be illustrated by adding a 
term — z/(V 2 ) 2 on the right of Equation (GQ), which results in a scaling of the roughness 
like 

w\L,t-Dr,v) = ^L^g>(^,^) . (9) 

Irrespective of t, the parameter tu/L 4 vanishes with increasing L much faster than 
the parameter tD/L 2 , so that w 2 (L, t; D, T 2 , v) approaches w 2 (L, t; D, T 2 , 0), assuming 
Q' \Dt j L 2 , tv j L 4 ) m Q(Dt/L 2 ) for sufficiently large L and irrespective of t, in particular 
irrespective of whether t is held fixed or the limit t — > oo is taken. This amounts to 
the statement that — z/(V 2 ) 2 is asymptotically irrelevant, so that the scaling of the 
Edwards- Wilkinson equation of the original form (TjQ) is recovered. 

2. The Edwards Wilkinson equation with drift 

Whether a term is deemed relevant is a matter of the canonical dimension of the coupling. 
Changing ([ID to the one- dimensional Edwards- Wilkinson equation with drift (EWd), 

d t (f>{x, t) = Dd?(f){x, t) + vd x <p(x, t) + rj(x, t) , (10) 

at first seems to suggest that D is irrelevant at any finite v, suggesting a scaling of the 
form 

w 2 (L jt] v,T 2 ) = ^- d g''^ , (11) 

so that 

a = (3 = z = l (suspected EWd) . (12) 

This result, however, is obviously flawed, if periodic boundary conditions (PBC) 
are applied. In that case the drift term can be transformed away by a Galilean 
transformation, <fi(x, t) = 4>(x—vt, t), so that <fi(x, t) follows the the original EW equation 
(jrj and thus the roughness of the interface displays the scaling derived in (JTj), 

a = 1/2 = 1/4 z = 2 (EWd with PBC) . (13) 

Since the initial guess (jl2"]) is based on a purely physical argument (rather than a 
mathematical one), this result merely questions the validity of this type of reasoning. It 
is clear that if v does not actually enter into the observable w 2 as defined in ([3]) because 
it can be transformed away, then it cannot dominate its scaling. As done below, where 
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v can not be transformed away, one might equally argue that D never becomes truly 
irrelevant, i.e. the exponents f fl3|) are a result of D being dangerously irrelevant [6]. 

A priori nothing can thus be said about the relevance of the couplings D and v. 
The scaling of the roughness has to be written as 

^i;ArV) = ^(|i) ( 14 ) 

which no longer fixes the scaling exponents a, (3 and z, as defined through Equation 
([6]), based on a scaling function on a single argument. The problem is the appearance 
of the new dimensionless parameter vt/L, which can alter the asymptotic behaviour of 
w 2 in a completely unknown way. In fact, with Dirichlet boundary conditions (BC), 
(p(x = 0,t) = 4>(x = L,t) = 0, it was found jS] that Q(Dt/L 2 ,vt/L) scales like 
(D I '(yL)) 1 / 2 as t — > oo at fixed L and like (Dt/ L 2 ) 1 ^ 2 for L — > oo at fixed t, in summary 

w 2 (L,t ] D^,v) = r 2 ^g>^ , (is) 

with Q\x) — > const, as x — > and Q'{x) oc y/x as x — > oo (with corrections in powers 
of Dt/L 2 ), so that 

a = 1/4 (3 = 1/4 z = l (EWd with Dirichlet BC) . (16) 

Depending on the boundary condition, the additional scale v thus either leaves the 
scaling of the Edwards- Wilkinson equation unchanged (as seen in the case of periodic 
BC), or changes them to anomalous values, which cannot be recovered from a simple 
dimensional analysis (as in the case of Dirichlet BC). The remainder of the present 
article is dedicated to the question which scaling behaviour is generated by Neumann 
boundary conditions, d x (f>(x = 0, t) = d x (p(x = L,t) = 0, which could lead, in principle, 
to a third set of exponents. 

2.1. Neumann boundary conditions 

If 0(x, t) is the height of a grown surface at time t atop a substrate at position x, the 
Dirichlet boundary condition mentioned above seems slightly unphysical, because it is 
difficult to imagine a mechanism that would pin the height to a certain value at the 
boundary of the surface, by taking up any excess amount of matter and providing it 
when needed. A zero flux boundary condition is much more naturally implemented, 
as it corresponds to imposing that no matter can enter or leave the substrate at the 
boundaries. The flux in (flO|) is 

j{x,t) = - {Dd x + v) <f>{x,t) (17) 

such that d t (f) = —d x j as a matter of mass balance (ignoring the noise). The condition 
j(x,t) = at x = 0, L is, however, very difficult to analyse. Since the true motivation 
of this study is the question to what extent the Edwards- Wilkinson equation with 
drift displays universal behaviour, in the following the weaker Neumann condition 
d x (j)(x, t) = for x = 0, L will be analysed. 
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To ease notation, it is helpful to absorb the various couplings and other dimensionful 
quantities in fflUj) into a redefinition of time, space, field and noise. With y = x/L G 
[0, 1], t = Dt/L 2 and coupling q = vL/D Equation ( flOj) becomes 



d r <p(y, r) = d y (p(y, r) + q<p(y, r) + £(y, r) 



where 



so that 



<p(v,r) = Jj^2<t>(x,t) and £(j/,t) = J j^n(x,t) (19) 



(ay,rny',r / )} = 5(y-y')5(r-T / ) . (20) 
The Neumann boundary conditions correspond to 

d y \ y=OA <p(y,T) = o. (21) 

As suggested by naive dimensional analysis, the coupling q seems to play an ever 
increasingly important role, as it diverges in the thermodynamic limit L — > oo. An 
alternative re-parameterisation of the drift that does not suffer from this problem is the 
dimensionless quantity u = tv 2 /D, which will be of great use below. 

The formal solution of Equation ffl8|) for a given realisation of the noise £(?/, r) 

( P(y, T )=[ dr' ( dy'G(y,T-r';y',q)^(y',r') (22) 
Jo Jo 

is based on the Green function G(y, r — r'; y', q) which describes the propagation of a 
Dirac delta peak at y' at time r' to y at time r. It is determined by considering the 
deterministic, homogeneous partial differential equation (PDE) 

d T G = (d 2 y + qd y )G (23) 

with initial condition lim^o G(y, r; y , q) = 6(y — yo) and Neumann boundary 
conditions, which in turn can be constructed from a complete set of eigenfunctions 
of the Sturm-Liouville problem 

K9n(y) = (d 2 y + qd y )g n (y) . (24) 

The operator can be made self-adjoint with the help of a suitable weight function, 
as discussed in great detail in [13], (f\g) = J^dy e qy f(y)g(y)- The set of normalised 
eigenfunctions orthogonal with respect to this scalar product is then found as 



9n(y) 



- qy . 



k n cos(k n y) + ^qsm(k n y) 



(25) 



V K + (V4)g 2 

and A n = — (l/4)g 2 — k\ for n > 1, where k n = irn. The eigenfunction without a node 
n = 0, deviates from that pattern, 



and Ao = 0. Having introduced a scalar product above that renders the differential 
operator self-adjoint, the temporal evolution of any initial distribution under the PDE 
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can be determined. On this basis, or equivalently, on the basis of completeness, 
the Green function is found to be [HI p. 63] 

G(y, r; y , q) = e™-?— + e^™)-^ £ * 



e« - 1 ti K + (V4)g 2 

x (k n cos(k n y ) + -q sm(k n y )^J (k n cos(k n y) + -q sin(& n y)^ (27) 

2.1.1. Poisson summation In principle, the Green function (1271) could now be used in 
(T22]) and the asymptotic properties of the roughness ([3]) be determined. However, it 
very quickly turns out that the real space limit of large L are very difficult to handle 
in Fourier space, k n , and the sums appearing in (1271) are practically intractable. A 
Poisson summation could be performed at any stage, but is most easily done directly 
on G(y, r; y , q) itself. To this end, one notes that 

f k n cos(k n y ) + \q sm(k n y )] (k n cos(k n y) + \q sin(k n y) 



kl + (l/4)q* y^—v^> ' 2 *~~vnwj ™ ■ ^ 

~(cos(A; n (2/ - y )) + cos(k n (y + y )j) 

+ \ k2 ~ (^4^2 sin ( k n(y + Vo)) - \q cos(k n (y + y )fj , (28) 

so that 

G(y, r; y Q , q) = e™-?— + e^^'K^ 
e q — 1 

oo 

x e ~ k " T ( cos ( k n{y ~ Vo)) + cos(k n (y + y )) 

n=l 



+ e -^q{y-yo)-jq 2 T 



qe 



n= 1 *S + (l/4)g 2 

x ^A; n sin(& n (y + y )) - -q cos(k n (y + j/ ))J • (29) 
This might not look like a significant improvement unless one notes that [151 P- 373] 



1 



oo 



3 Cz, r) = l + 2^ e mTn cos(2n7r^) = -_ ^ e" 1 -^ (30) 

n=l y—1T n= _oo 

which allows the re-summation of the two cosine terms in ( 129]) . 

x (k n sm(k n (y + y )) - ^ cos(fc n (y + y ))) + e-^-*°H^-J= 

^ / (y-j/Q + 2n) 2 fa + H0 + 2n) 2 \ 

x 22 (e ~r +e J? J . (31) 

n=— 00 ^ ' 

Without the exponential prefactor, exp[— (l/2)q(y — y ) — (l/4)g 2 r], the last summation 
describes the diffusive motion of particles with (positive) mirror image on a ring of 
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circumference 2 not being subject to the drift. In the following, the aim is to express 
Equation (l3Tj) in terms of the free propagator on R 

1 (v + Qt) 2 

$ ( 2/ ,r;g) = -=e-T , (32) 

which solves f[2"3"j) . so that none of the terms is expressed any longer in Fourier space, 
which facilitates integration and the determination of the asymptote using a saddle point 
approximation. In fact, including all pre-factors, the last sum in ( 13 ip can be written as 

oo 

J2 e^($ (y - Vo + 2n, r; q) + e™$ (y + y + 2n, r; q)) . (33) 

n=— oo 

The first three terms in (1311) . the two exponentials and the sum, are much more 
difficult to handle. Taking the limit r — > 0, the expression 



s ^ yo} ~ e #-i +e hK + iVW 

x (k n sm(k n (y + y )) - -q cos(k n (y + y ) )J (34) 

ae qyo 1 i , % 00 1 

2 g A£ + (1/4)9* 

x ^A; n sin(& n (y + y )) - cos ( k n(y + Vo))) (35) 

can be regarded as an initial source which is propagated under the PDE (T23|) and on R, 
by the propagator ( l32l) . so that in fact 



G(y, r; y , ?) = / dy s(j/, 2/0)^0(2/ - & r; g) 

00 

+ Yl e 9n ($ (y - yo + 2n, r; g) + e «°$ (y + y + 2n, r; g)) . (36) 

n=— 00 

Taking the limit r — >■ for the entire G(y, r; y , q), (15111 reveals 

00 

hmG^, r; y , g) = s(y, y ) + e -3*(v-w) £ (% ~ ^ + 2n) + % + y + 2n)) (37) 

n=— oo 

which means that s(y, ?/ ) = for y, y e [0, 1], since lim^o G(y, r; y , q) = 5{y — yo) on 
that interval, which implies that 

00 1 / 1 \ 2 

(1/4)^ r- sin(M - 2 gcos(M ) = -^i e ^ 2 (38) 

for any z G [0, 2]. Clearly, the left hand side is periodic in z with period 2, which means 
that it is in fact the periodic repetition of the right hand side, i.e. 
00 1 / 1 ^00 

^ ui/i iA\ n 2 (knsm{k n z) - -qcos{k n z)j = - £ p(z + 2n; q) (39) 

n=-oo K n 1" I 1 / 4 / 1 ? V Z / n=-oo 

with 

/)(z;?)~e^ 2 / M (z) (40) 
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and the indicator function 

T , . f 1 for z E ft , ... 

Ia(z) = { otherwise . (41) 

On this basis, the source can be rewritten as 

*(y.W)) = ^3j;(l- E e^J 0j2 ( 2 / + |/o + 2n)J (42) 

and therefore 

/•2 (/e Wo 00 

G(y,T;y ,q)= / # J ^ T E ~ e 9 ")$o(l/ + J/o + 2n - y, r; q) 

Jo e l n= _ 00 

oo 

+ E e qn (® (y-yo + 2n,T;q) + e< i y o <S>o(y + yo + 2n,T;q)) . (43) 

n=— oo 

Problems of convergence have not been addressed here in a detail, which could affect 
some of the manipulations done above, in particular when the order of integration and 
summation is changed. Yet, because of the free propagator $o effectively posing an 
exponential cutoff on the sums as well as the integrals, all convergence issues turn out 
to be harmless. 

To verify that the propagator ( l4"3j) indeed solves the PDE ( 1231) is a matter of 
a straight forward calculation. Similarly, the initial condition lim^o G(y, r; y , q) = 
S(y — 7/o ) f° r y-, Vp ^ [0, 1] can be identified by mere inspection. The only slight difficulty 
are the boundary conditions d y \ y= o tl G(y, r; y , q) = 0, which are most easily verified by 
writing the gradient as 



d y G(y,r;y ,q) = e ^ y Vo) 4 qT -==i ^ e *r - e 

V47TT I 2r 



' e ^ + e & > (44) 



2r 



Calculation of the roughness 



The propagator (143!) can be used in the formal solution (1221) . which provides the basis 
for calculating the roughness via ([2D and (TP9"j) . 

w 2 (L, t) = — — / dyi / dy 2 [t%i - t/ 2 ) - r)^(y 2 , r) , (45) 



D Jo * Jo 

where the ensemble average enters only through ( 1201 ). 

The integral on the right is the dimensionless roughness, 
i ,i 

d7/i / dy 2 [5(t/i - 1/2) - Imyi, T)(p(y 2 , r) , (46) 

which depends on only two dimensionless parameters, r = Dt/L 2 and g = vL/D. 
However, as discussed in the following, the parameterisation 

w 2 (q,u) = q[ dyi f dy 2 [S(yi - y 2 ) - l](p{yi, r)p(y 2 , r) , (47) 
Jo Jo 
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with u = tv 2 /D = q 2 r is advantageous, so that w 2 (L,t) = (T 2 /v)w 2 (q,u), which means 
that all scaling of the roughness can be read off the scaling of w 2 without any further 
prefactor such as L in f j45l) . 

The limit t — > oo in ( 145]) translates simply to r — > oo in ( H6|) and the scaling of this 
limit in L is observed through the scaling in q, which is unaffected by the limit r — > oo. 
So, in order to determine the roughness exponent a, the parameterisation of ( 146]) in q, r 
is suitable. However, taking the thermodynamic limit L — > oo means r — > and q — > oo 
simultaneously and not independently, since q 2 r = u remains constant. 

Using w 2 (q = vL/D,u = tv 2 /D) as the parameterisation of the roughness means 
that the two limits L — > oo and t — >• oo affect each only one dimensionless parameter 
but not the other. The limit t — > oo means u — > oo and the scaling in large L is 
observed through arbitrarily large but fixed q. The limit L — > oo means q — > oo and the 
scaling in large t is observed through the scaling in large but fixed u. The two different 
asymptotes will be determined by considering large q in both cases and u/q diverging 
or u/q vanishing. 

Using the parameter u in favour of r needs to be implemented in the propagator as 
well, ( H3|) . using r = u/q 2 . It is also useful to replace the integral over r in the formal 
solution (|22|) correspondingly by an integral over u so that 

v 

— w 2 (L,t) = w 2 (q,u) = q- 1 J dy 1 dy 2 \5(y 1 -y 2 ) - 1] 



x 



/ dy' / du'G(yi,u';y',q)G(y 2 ,u';y',q) , (4£ 

JO JO 



where G is the re-parameterised propagator with all parameters explicitly appearing as 
arguments. 

The bulk of the work is performing the integration ( HBl) . This can be done more 
conveniently by splitting the propagator up into four terms, 

G{ yi ,u';y',q) = A, - B, + Q + D, , (49) 



where 



qe qy 



A,, = -f— - (50a) 
e q — 1 



B, = - £ / dy* $0(1/, + y' + 2n t - y h u; q) (50b) 

00 

Q = £ e^^o^ - y' + 2th, u; q) (50c) 



with 



Di= £ e^ + ^$ (y l + y , + 2n,, M ;g) (50d) 

rij=— 00 

(50e) 



$o(y,w;g) = -?=e *» . (51) 

V47TU 
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2.3. Integration and analysis 

Because of the prefactor [S(yi — y 2 ) — 1], for each of the ten distinct terms generated by 
using ( l50i) in f HHi) . effectively two different integrations have to be performed, namely 
one over yi and y2 distinct, and one over y\ = y 2 . In addition, an integral over y' and 
u' needs to be performed, as well as over y for term Bj. In all saddle point 

approximation (SPA) is used, specifically 

-f= dye « = J A h/o - - for u/q 2 < 1 , 52 
V 4vtm M \ q J 

where Ia{u) is again an indicator function, (|4ip . It is worth stressing that there are no 
further algebraic terms, i. e. all other contributions from the integral vanish exponentially 
in large q 2 /u. 

In most cases, further integrals over variables contained in the argument of the 
indicator function are to be performed. If the intersection of the integration range of 
the other variables and the set A above has finite measure, the result is straight forward 
to calculate. If, however, the intersection contains a single point, then the SPA has to 
be considered to higher order. In general, for each "marginal variable" the resulting 
integral acquires an additional factor y/u/q. For example, to leading order u/q 2 , 



f 1 i f 1 , Mvi-v2-j/3)) 2 \/4nu f 1 , T . 1 . . Ayfnu , . 

/ dyi / dy 2 e w / dy 2 Vi) W + 1 / 3 ) = ~^ — ( 53 ) 

Jo Jo q Jo oq 

using ( |52j) in a non- marginal case (maximum of the exponential at y± — y 2 = 1/3 G 

(—1,1)), whereas replacing 1/3 by 1 produces 

i /- 1 , (g(«i-» 2 -i)) 2 2m . a . 

dyi dy 2 e ^ w — , (54) 
o Jo g z 

rather than ;y ^^ J 1 d7/ 2 1(o,i)(z/i + 1) = 0. Integrals like that result in subleading terms, 
whose amplitude is not normally calculated in the following. Instead, only the power of 
q is noted and by comparison to other terms it is verified that it is safe to ignore it. The 
same result is recovered by power counting in the four terms fl50l) . where each integral 
gives rise to a leading order q' 1 and thus each of the four terms has the same algebraic 
dependence on q. 

Further details of the calculation are exemplified in the appendix. Combining all 
contributions gives 

T 2 r 1 r 1 
= — Q' 1 / dyi / dy 2 [<%x - y 2 ) - 1] 
v Jo Jo 

x (Ai - B x + Ci + Di) (A 2 - B 2 + C 2 + D 2 ) 
For u small, looking at the limit q — > oo, to leading order: 

lim w 2 , . , = — \f— (55) 

(J^-OO Physical v V 2lT 

For q large and considering the limit u — > oo, to leading order: 

T 2 2 

lim w 2 . . , = i=\[q (56) 
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3. Discussion and conclusion 



Upon replacing u and q by their definitions, u = tv 2 /D and q = vL/D, the two 
asymptotes for the roughness derived in Equation fl56|) and fl55|) are 



tt-^h for L — > oo 

w*(t,L) = T 2 x { (57) 
I\/r4r for t -> oo 

3 V 2nDv 

which is, to leading order, identical to the result in [5] for Dirichlet boundary conditions. 
The exponents as defined in ([6]), 

a = -, $ = z = 1 (EWd with Neumann BC) (58) 

therefore reproduce (fl6l) . On the one hand, this is a surprising result, because after 
realising that the usual dimensional arguments do no longer apply, any exponents 
are mathematically possible. The fact that Neumann boundary conditions reproduce 
the anomalous results of the Dirichlet case (contrasting those for periodic boundary 
conditions) even down to the amplitudes therefore point to some universal mechanism 
at work in both equations. Physically, on the other hand, arguments very similar to 
those discussed in [5] apply: The drift with velocity v effectively constantly re-initialises 
the interface with vanishing slope from one side to the other, while constantly under 
the influence of the external noise. This mechanism erases practically all features that 
develop over times exceeding L/v and thus in saturation, t — > oo reduces the roughness 
to about the value after time L/v. In fact, the stationary roughness is 2/3 of what is 
extrapolated from initial roughening: 

lim w 2 (t,L) = \ lim w 2 (L/v,L') (59) 

In conclusion we have shown that the Edwards- Wilkinson equation with drift 
and Neumann boundary conditions produces anomalous scaling different from what 
is expected from naive dimensional analysis and easily derived for periodic boundary 
conditions. The scaling and the amplitudes of the Dirichlet case are reproduced, 
consistent with a simple physical scenario of an interface that is constantly re-initialised. 
Neumann boundary conditions are physically more relevant than Dirichlet ones, yet far 
more difficult to handle analytically. 
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Appendix A. Integral ( I48j) term by term 



In this appendix some details of the integral ( 1481) are shown, considering the propagator 
term by term as defined in (j4"9j) and fl50l) . 

Appendix A. 1. AxX 2 Term 

Since there is no y± or y 2 dependence in Ai, (150aj) : 

/>1 pi pi pu 

/ dyx / dy 2 [5( yi - y 2 ) - 1] / dy' / du' A,X 2 = 
jo Jo Jo Jo 

=0 

for X being any of the A, B, C or D terms. 
Appendix A. 2. BxXj Term 

Since all BxXj integrals follow similar arguments, we exemplify the procedure on BiC 2 . 
The integrals to be calculated (to leading order) are 

<c = —Q' 1 /\i /W2 [%i - y 2 ) - 1] / V f du' (A.l) 

f JO JO JO JO 

x — \ e g 1 / dyi . e > e q 2 . e *« 

We consider the total exponential of B1C2 in the form 

exp(q(m + n 2 ) - f(y',y)) (A.2) 



e« - 1 
where 

f(y',y) = -qy' + ^L+y' + 2 ni -y + ^j +^(y 2 - y > + 2n 2 + ^J . (A.3) 

A contribution in the form (1A.2|) vanishes for large q unless 

n x + n 2 - IMllL > 1 . (A.4) 
Q 

Applying the SPA to the integral over y', means to evaluate the integral at y' that 
minimises /: 

1 1 u' 

y'o = -(j/a -Vi) + ("2 -nx) + -y + — , (A.5) 



where / now becomes: 



f(y'o, y) = ^ + 2/2) + (nx + n 2 ) - I) + g(l/i + 2rn - y) . (A.6) 

Applying the SPA again to the integral over y gives, with f(y' ,y) above, 

Au' 

yo = yi + y 2 + 2( ni + n 2 ) + — , (A.7) 

Q 
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where now 

f{y' ,y ) = -q{y 2 + 2n 2 )-2u' . (A.8) 

Determining the asymptotic behaviour is greatly facilitated by writing — = a + m > 0, 
with a G [0, 1) and m G IN , so that \u' /q\ = m. Using the final result of /, (lA.8j) . the 
inequality in ( ]A.4j) becomes: 

f 

ni + n 2 - - = y 2 + 2a + m + 3n 2 + 2m > 1 . (A. 9) 

Q 

As y 2 G [0, 1] and a G [0, 1), so y 2 + 2a G [0, 3), the contribution (IA.20 vanishes for large 
q unless 

n\ + 3n 2 + 2m > —2 . 
Since ni, n 2 and m are integers, we have the inequality 

ni + 3n 2 + 2m > -1 . (A. 10) 

Both SPAs produce an indicator function as well. Writing — in terms of a and m gives 

I(o,i) (l/o) = I(o,i) [(ya " yi)/2 + (n 2 - th) + y/2 + a + m] (A.ll) 
I(o,2) (yo) = I(o,2) [yi + 1/2 + 2(ni + n 2 ) + Aa + 4m] . (A. 12) 

From the indicator function of y' above, (IA.11I) . we have 

< ^(y 2 -yi) + (n 2 -n 1 ) + ^y + a + m<l (A.13) 

as a condition for ( 1A.2I) to contribute. As yi,y 2 G [0, 1), y G [0,2) and a G [0, 1), this 
implies that \{y 2 — yi) + \y + o G [— |, |). By subtracting this range, (1A.13|) becomes 
— | < n 2 — ni + m < |j which gives 

-2<n 2 -n! + m<l (A.14) 

as ni,n 2 and m are integers. Similarly, for flA.12j) to be non-zero, we have the condition 

-2 < (m + n 2 ) + 2m < . (A.15) 

BiC 2 vanishes if either one of these two ranges, f ]A.14j) and (1A.15|) . is incompatible 
with all possible values of ni,n 2 and m. By adding ( lA.lOj) and (1A. 14j) we find 
-3 < 4n 2 + 3m and as m 6 H° it follows n 2 > 0. Adding (IA.14|) and (1A.15[) . we 
obtain —4 < 2n 2 + 3m < 1 and with m > it now follows n 2 = m = to prevent either 
indicator function from vanishing. Therefore, the two indicator functions become 

1(0,1) (y' ) = 1(0,1)13/2 + 3(T ] ( A - 16 ) 
1(0,2) (y ) = I(o,2) [yi + y2 + 2n, + 4a] . (A. 17) 



The next variable to determine is rii. With n 2 = m = Equation flA.9j) reads: 

y 2 + 2a + n 1 >l (A.18) 

Multiplying flA.18[) by 2 and rearranging, we have the condition y 2 +Aa+2n x > 2—y 2 > 1, 
as y 2 G [0, 1), and thus 

y 2 + 3a> \-o-2n x (A.19) 
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Using n 2 = m = in ( lA.lOj) and HA.15[) gives the range of n! as — 1 < < 0, whereas 
(IA.16I) indicates 

1 > y 2 + 3(7 > . (A.20) 

If rii = — 1, then (1A.19I) becomes y 2 +3a > 3 — a. Since cr g [0, 1), this contradicts (1A.20j) . 
Therefore, we conclude that n\ = 0. However, if we use the fact that n\ = n 2 = m = 
in (1A.9I) we get y 2 + 2a > 1. By comparing to ( 1A.20I) . we have a < which contradicts 
the fact that a G [0, 1). As the conditions cannot be fulfilled simultaneously, this implies 
that the integrals vanish exponentially in large q. 

Considering in addition any marginal cases will produce terms of lower algebraic 
order in q. For the present integral, the marginal cases are n\ = n 2 = m = with 
y 2 = and a = (double marginal), as well as n\ — 1, n 2 — m — with y\ — y 2 — 
and cr = (triple marginal). Power counting in the initial integral flA.lj) thus gives 



overall contributions to w% c of order 0(q~ 2 ) and 0(q~ 3 ) respectively. Extra care must 
be taken when considering the integral with upper bound u, as u — > oo might be taken 
before q — > oo, in which case the integral over u' might give rise to a term of order q 
itself (see, for example, flA.261) versus flA.271) ). In the present case this does not apply, 
because u' = q{m + a) and both m as well as a are fixed by the SPA. 
By similar arguments, the terms B1B2 and B1D2 vanish. 

Appendix A. 3. C{K 2 and DiD 2 Term 

Calculations for the terms C1X2 and D1D2 are very similar and we therefore exemplify 
the procedure for D!D 2 only. 

We write contribution by D!D 2 as 

r 2 

o J- 



^DD 



-q- 1 / d Vl / dy 2 [5{ yi - y 2 ) - 1] / dy' / du' 
v Jo Jo Jo Jo 

X q -lq<,Vl+v' ' + 2n 1 ) + u'] 2 . ^ O - [gfag W +2n 2 ) W] 2 



x V e gni ^=e w e ^ V e gn2 ^=e ' ' w J e TO ' 

We consider the total exponent of D!D 2 in the form 

exp(g(m + n 2 ) - /) (A.21) 
We first apply the SPA to the integral over y', and obtain the minimum y' , so that 



2 



^(yi -y-i) + (m -n 2 )} 



2 



+ 2q 



^(2/1 + 2/2) + (ni + n 2 ) 



Again, the analysis of the asymptotic behaviour is greatly facilitated by writing 
— = a + m, with o G [0, 1) and m G IN , so that the indicator function for y' Q is 

l (o,i)(y'o) = I(o,i) (°" - \{yi + I/2) +m- (m +n 2 )^) . (A.22) 

Since cr — |(|/i + 2/2) £ [ — 1> 1), the remainder of the argument m — (ni + n 2 ) must be 
greater than —1 or less than 2 for this term to contribute. Therefore we will only need 
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to consider only n\ + n 2 = m and n\ + n 2 — m — 1, incorporated below by means of 
Kronecker 5-functions. The total contribution then becomes: 

^2 



^DD 



r r 1 r 1 
—q' 1 / dyi / dy 2 [<%i - y 2 ) - 1] 
v Jo Jo 



X 



ru 

L M s 

JyJ ni,n,2 



V87TU 



= exp 



2u' 



n 2 



-Xyi-yi) + (ni -n 2 ) 



-2q 



;(yi + 2/2) + (ni +n 2 ) 



e- qm ko,i) 



} m,n\+ri2 



+e - ? (l-m)j ( _ i| 



0) 



J m- l,ni+n2 



At this stage, it is sensible to consider the 5-term and 1-term in [5(yi— y 2 ) — 1] separately. 
For 2/1 = 2/2= 2/j i- e - the contribution of the 5-term, the integral to consider is 

hi I dn! — ran I — 

2u' 



(DxD 2 )i = J\ Y, ^=exp(-|l(n 1 -n 2 ) 2 -2g(2/ + (n 1 +n 2 )) 



711,712 



X 



e 9m I( ,i)[o- - 2/]5 m ,„ 1+ „ 2 + e 9(1 m) I(-i,o)] [<? - y]5 m _i >ni+n2 

using I(o,i) (y + 1) — I(-i,o)(y)- This term can in turn be split into three parts, 
each accounting for one of the Kronecker (^-functions and the possible values of m. 
To prevent exponential suppression, m = is required for the prefactor e~ qm and 
m — 0, 1 for the prefactor e~ q ^~ m \ which gives rise to three terms, (D 1 D 2 )«5 = 
(D 1 D 2 ) (5l + (D 1 D 2 ) ( 5 2 + (D 1 D 2 ) ( 53. 

In the first term (D 1 D 2 )s 1 we have m = and n\ = —n 2 from the Kronecker 5. 
Ensuring the indicator function is non-zero, we change the integration limits of y and 
using 5o,ni+n 2 we write n\ = —n 2 = n, 

(DiD 2 )<5 = fdy f U du'Y^=^(- q ^ 2 )-2qy] 
Jo Jo „ v8vtm' V u J 



q 



-2qa 



W ) V87TM 7 

In the limit of large q, the first exponential will asymptotically vanish unless n = 0. 
Since \u'/q\ = m = 0, we change the upper integration limit of u' to min(g, u), and 
apply n = 0: 

/■min(l,2) 1 / „2 \ i 

(DiD 2 ), 1= jf 9 ,da£^exp(-^(2n 2 ^ 



n 2g 



1 /-min(l,|) 




2tt Jo 



da 



2 min 1 , 



u 



1 -e 
1 



2qa 



1 - e~ 2qcT 



1 



r -, 2gmin 1, 



using -u' = ag and 

[ y dx x^e~ x = + 1,2/) , 
Jo 

which converges to ^/n for ji = 1/2 in the limit of large y. 
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Similar arguments can be used for (D 1 D 2 ),5 2 . In that case, the support of the 
indicator function changes from [0, 1] to [—1,0] since m = 1 and therefore the integral 
changes to: 

(D 1 D 2 ) S2 = [°°dy f U du' £ -L= i exp (-^{2n 2 ) - 2qy] 
Jo Jo n V8iru' V u J 



i f± 

4V27T 



=T [ - , 2omin f 1, — 
2q I 2' I q, 



Finally, the contribution of (D 1 D 2 )«5 3 , with m = 0, is 

(DiD 2 ), 3 = fdy f U du' £ -7= 
Jo Jo nun2 V8nu' 

x exp ("Ju 7 ^ 1 _ n2 ^ ~ 2q ^ y + ( Ul + n2 ^) e ~ 9l (- 1 ^ a ~~ v) ■ 

Since m = 0, this means that n x + n 2 = —1. Letting n\ = n, this gives n\ — n 2 = 2n + 1 
and thus 

(D,D 2)J , = J\ /V ? -^exp (-, (X(2„ + 1)» + 2y - 1 

xe^I(-i, )((T-!/) • • 
Each exponential in the sum vanishes for large q unless 



^(2n + l) 2 + 2j/-l<0. 



(A.23) 



Since m = 0, u' = qa < q, and from the indicator function follows y > a, so that 

> ^-(2n + l) 2 + 2y-l>^- + 2y-l>^- + 2a-l (A.24) 
2a 2o~ 2a 

which implies —2a > (2a — l) 2 for a > 0, which is impossible. The (DiD 2 )<5 3 term 
therefore vanishes. 

In summary we have 



[D 1 D 2 ) 5 = (D 1 D 2 ) Sl + (D 1 D 2 ) S2 + (D 1 D 2 )s 3 



2V27T 



min 1, 



u 



For u small and by taking the limit q — > oo, it is obvious that: 
lim(D 1 D 2 ) 5 ^0 

Yet, taking u — > oo first gives min (l, |) = 1 and we find to leading order in q: 

r 2 _ x i rg-_ r 2 1 



lim^D,)^-^ 1 -,/-^ = 



u^oo- ' v 2 V 2n v 2y/2nq 

Next we consider the roughness contribution of the 1-term, i.e. yi ^ y 2 . 



(A.25) 
(A.26) 

(A.27) 



:D 1 D 2 ) 1 = l dm I dy 2 [du> ^ 7= ( - ^ 



11,12 



iru' 



+e -9(m-l) I 



(-1.0) 



"2(2/1 + 2/2) +0" 



2w 

e" 9m I(o,i) 

^m— l,ni + «2 



2 rl 

2(2/1 -2/2) + (ni -n 2 )] 



-^(2/1 + 2/2) + o- 



+n 2 
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This term, again, can be split into three parts, (DiD 2 )i = (D 1 D 2 )i 1 + (D!D 2 )i 2 + 
(D 1 D 2 )i 3 . Using similar as for the 5-term, (D 1 D 2 )i 1 contribution vanishes for large q 
unless m — 0, which implies n\ = —n 2 = n. An SPA applied to the integral over yi 
gives a minimum at 

Au' A 
Vio = V2 - — - 4n , 



so that the exponent is 



2u! 



f(y lo ) = q(2y 2 - — -An). 



in the notation introduced above. The first term therefore gives 

(D 1 D 2 ) ll = [ 1 dy 2 £ f U du' 
Jo „ Jo 



— — exp \-q\2y 2 An 

q 2 /8u' V V q , 



xl 



(0,1) 



- 1 ( Au ' A \ 

-y I V2 ~ — ~ An + y 2 J + a 



(o,i) 



Au' A 
An 

Q 



rl ru 

= dy 2 du exp (-q(2y 2 - 2a - An)) 
Jo Jo 

x I(o,i) [-V2 + 3o- + 2n]I( 0) i) [y 2 - 4o- - An] 
Both indicator functions suggest n = or n = — 1, otherwise the (DiD 2 )i 1 term will 
not contribute. In the limit of large q, the term will be exponentially suppressed for 
n = — 1 so that the only case to be considered is n = 1. However, using n = in both 
indicator functions implies o < for this term to contribute which contradicts a > 0. 
We conclude that (DiD 2 )i 1 does not contribute. 



By observation, (D!D 2 ) l2 is identical to (D 1 D 2 ) ll with the domain of the indicator 
function changing from [0, 1] to [—1, 0] as m = 1. 

Q 



rl ru 

(D!D 2 ) l2 = / dy 2 Yl / d«' 
Jo „ Jo 



7T 



7 exp 



x !(-i,o) 



-1 / 4u' \ 

2~ \ 2 ? An + y 2 j + a 



-q \ 2y 2 An 



(o,i) 



Au' A 
y 2 An 



dy 2 / du' exp {—q{2y 2 — 2m — 2a — An)) 
io Jo 

x I(_i )0 )[-y2 + 3(7 + 3m + 2n]I( ,i)[y 2 - 4<7 - 4m - An] 

rl ru 

= dy 2 du' exp {-q(2y 2 -2 -2a- An)) 
Jo Jo 

x I(_4,-3)[-2/2 + 3(7 + 2n]I( 4>5 )[j/2 - 4(7 - 4ra] . 
Both indicator functions contribute only if n = —2, in which case, however, the term 
vanishes exponentially in large q. 

Finally, we consider the (DiD 2 )i 3 contribution, for which m = and n\ + n 2 = — 1. 
By writing n = n\ = — 1 — n 2 and applying SPA to the integral over y\ we find the 
minimum of the exponent at 

Au! 



yi = V2 



2(2n + l)=y 2 -Aa- 2(2n + 1) 
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2u' 

/(yio) = ?(2y 2 - — - 4n - 1) = g(2y 2 - 2<r - 4n - 1) . 



and thus 
(DiD 2 ) l3 



/ d</ 2 £ / dvf 
Jo ~ Jo 



71 



TXU 



q 2 /8u 



-exp(-q(2y 2 - 2a - An - 1)) 



x I(-i,o) 



-1 



(y 2 - 4a - 2(2n + 1) + y 2 ) + a 



I ( o,i)[!/2 -4a -2(2n+l)] 



dy 2 / du exp {—q{2y 2 — 2a — An — 1)) 
•A) 

xI(-2,-i) [-2/2 + 3a + 2n]I( 2)3 ) [vi - Aa - An] . 

Similar arguments as for (DiD 2 )i 2 apply. Both indicator functions contribute only if 
n = —1, where the term is exponentially suppressed. Therefore, (DxD 2 )i 3 does not 
contribute either. 

Hence in total, considering both limits we have to leading order: 

r 2 i 



lim Wq D = and lim u> 2 



q— >oo 



DD 



v 2^2^q~ 



(A.28) 



Any marginal cases are necessarily subleading compared to (1A.25j) and thus can be safely 
ignored. 

By similar procedures, we obtain the leading order behaviour of Wq C and Wq D as 
follows: 



q— >oo 



lim w cc = — \l — and lim w 



v V 2tt 



cc 



V 3V27T 



(A.29) 



and 



lim Wen = and lim = . 



In summary, the only non-zero contributions are f ]A.28j) and ( 1A.29|) . 
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